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Abstract 

The convergence characteristics of two viscous core corrections as used 
in straight-line segmentation methods are rigorously analysed. These are 
curvature corrections that account for the induced velocity contribution 
at a point on a vortex filament due to local curvature and core corrections 
that remove unrealistically high induced velocities near vortex segments. 
Two alternative versions of the latter correction are studied: the original 
as introduced by Scully and a recently introduced improved correction by 
the authors. The problem is analysed using a viscous vortex ring. A high- 
order numerical filament method is presented that uses rational B-spline 
curves to model the vortex ring geometry exactly. First, the separate in- 
fluence of the two corrections on induced velocity values are studied using 
the numerical filament method. Afterwards, the corrections are combined 
in four separate cases that are studied using the segmentation method. 
It is found that in order to get accurate results at coarse discretisations, 
curvature corrections cannot be neglected. When used together with the 
original core correction model, results start to diverge from the reference 
values for discretisations smaller than approximately ten degrees. Com- 
bining the curvature correction with the improved core correction extends 
the discretisation region for which accurate results are found down by ap- 
proximately one order. Finally, it is shown that segmentation methods 
that use the original core correction model without curvature correction 
completely fail to converge to the reference values. 



* Submitted for review to the Journal of the Ameriean Helieopter Society on November 23, 
2011. 
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1 Introduction 



Engineering rotorcraft simulations that include a free wake component mainly 
use straight vortex segments to discretise the curved vortices that are trailed 
from the rotor blades (Refs. [ll|, 2^ 2^)- An analytical solution of the Biot- 
Savart law is available for this basic wake element (Ref. which is one reason 
for its frequent use. A disadvantage of this method is that the analytical solution 
must be corrected to remove unrealistically high induced velocity values close 
to the segment. 

A consequence of the above approach is that for a marker point on a curved 
filament, the contributions due to the two directly adjacent segments are zero. 
In literature, there is disagreement over the need to correct for this contribution 
at marker points. Bagai (Ref. fS]) and Bhagwat and Leishman (Ref. [^) claim 
it is not necessary, while others (Refs. [14, 1^ 2^) claim or implicitely assume it 
is a contribution that should not be neglected. Other researchers acknowledge 
that both core corrections and curvature corrections are needed, but then only 
use the former (Ref. [9]). 

The core correction model introduced by Scully (Ref. is used in rotor- 
craft wake simulation software to remove unrealistically high velocities near a 
filament. Recently (Ref. [IHl), it was shown that when curvature corrections are 
neglected in induced velocity computations (as in Refs. 0,3), the use of Scully's 
core correction model results in an underprediction of the induced velocity of a 
viscous vortex ring by approximately 40%. An improved core correction model 
was presented that reduces the error considerably at fine discretisations for a 
large range of relative core sizes. The reference formulas used in Ref. 26[ were 
those for a vortex ring with a non-circular core section (Ref. [l^]), while in nu- 
merical filament methods, a circular core shape is assumed. Recently (Ref. [25|V 
it has come to our attention that this assumption may lead to different con- 
stants in the analytical formulas for the induced velocity of viscous vortex rings 
(Ref. 0,113). 

The conclusion in Ref. Q that curvature corrections are superfiuous at equiv- 
alent arc lengths smaller than approximately 10 degrees is based on the use of 
the original core correction model. A similar conclusion may not hold when the 
improved correction is used. 

Therefore, the purpose of this paper is to investigate the convergence char- 
acteristics of straight-line segmentation methods when used with 1) the original 
core correction without curvature correction, 2) the original core correction with 
curvature correction, 3) the improved core correction without curvature correc- 
tion and 4) the improved core correction with curvature correction. Using this 
investigation, we would like to give a definite answer to the question whether or 
not curvature corrections are necessary in straight-line segmentation methods. 
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2 Induced Velocity Computations using the Biot- 
Savart Law 



2.1 Desingularisation of the Biot-Savart Law with Veloc- 
ity Smoothing Functions 

The velocity at a point P induced by a singular vortex filament C(u) with 
strength T is given by Leonard (Ref. ^ISj) as 

r f rp — r{u) ^ dr{u) ^ 

^ 47r Jc(u) |rp - du 

= T [ K{vp-v{u))-K^^du, (1) 
Jc(u) du 

Av{u)du, 

C{u) 

where is the derivative at a parametric point u on the filament with respect 
to the cm^ve parameterisation, K(r) — (— l/(47r|rp))r is the Biot-Savart kernel, 
r = rp — r(u) and Tp is the position vector of the point P. 

In the vortex dynamics literature (Refs. 0, H^l), two methods are used 
to desingularise the Biot-Savart law, the cut-off method and the smoothing 
method. In the former, a portion of the integral near the point is removed 
from the integration domain, which removes the singular part. The width of 
the region that must be excluded is proportional to the relative core size a / R. 
The exact width depends on the vorticity distribution in the core. For example, 
for a core with a uniform distribution of vorticity, the proportionality constant 
can be found to be ^e^/'* (Ref. 22], p. 213). This method can only be used 
for (parametric) points located on the filament and will not be used here. In 
the smoothing method, the singular Biot-Savart kernel K(r) in Eq.[T]is replaced 
with a mollified kernel 'K.i,{r) which removes the singular behaviour at the vortex 
center line, 

K^W--^-' (2) 

where giy{v) = g {^-^^ — g{p^ is a three-dimensional velocity smoothing func- 
tion (Ref. [2^) with cr the core size and p the relative distance of a point with 
respect to a point on a filament. The regularised Biot-Savart law can now be 
written as 

vp^r/ K.(rp-r(.))x^rf..-r/ x ^d., (3) 

Jc(u) du Jc(u) du 

which can be integrated numerically using any suitable quadrature method. The 
method used to describe the filament geometry is discussed later. 

Related to the velocity smoothing is the vorticity smoothing function Ccr that 
yields the vorticity field (Ref. Q) 

c^.(r) =T [ Urp - r{u))^^du. (4) 
JC{u) c/u 



^Whon a smoothing function is displayed without a subscript to denote the dimensionality, 
it is assumed three-dimensional. 
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where Co-lr) — ~ ^C(p)- Given a three-dimensional vorticity smooth- 

ing function C3(p)i ttis associated three-dimensional velocity smoothing function 
53 (p) can be found from (Ref. [29l] ) 

9z{p)= \\z{t)edt. (5) 
Jo 

For an infinitely long, straight filament, the induced velocity field is circum- 
ferential. The two-dimensional velocity smoothing 172 (p) that induces the same 
velocity field is obtained from (Ref. [28|). 

g,{p) = r ^ dt. (6) 

Similarly, the vorticity distribution in two dimensions Qii^p) can be found from 
(Ref. (3), 

In two dimensions, the velocity and vorticity smoothing functions are related 
through g2(p) = Jq C2{t)tdt. Lastly, the normalised swirl velocity profile vg{p) 
can be found from the two-dimensional velocity smoothing function as followqj, 

. ^ 27r.g2(p) 

Mp) = • (8) 

P 

It is apparent from Eqs.[n]and[7]that computing the two-dimensional smoothings 
from the three-dimensional smoothings is straightforward, but the opposite can 
be very difficult. 



2.2 Example Velocity Smoothing Functions 

Winckelmans and Leonard (Ref. [29] ) give an extensive list of two- and three- 
dimensional vorticity and velocity smoothing functions with their convergence 
characteristics. Nowadays, these smoothings are mainly used in applications 
where a very large set of vortex particles is used to model regions of vorticity. 
The exact swirl profile induced by a single regularised particle is not as impor- 
tant as the vorticity distribution of a large set of particles. However, when a 
single parametric curve is used to model a single curved vortex filament (as in 
helicopter applications), one specific characteristic of the swirl velocity profile 
is important, namely the location of the extremum. As an example, the swirl 
velocity profile associated with the high-order algebraic smoothing introduced 
by Winckelmans and Leonard (Ref. [29|) does not have its extremum at the 
nondimensional distance p = 1, but more inward (Ref. j25j'). For this reason, it 
is not used here. 

Some smoothings relevant to rotorcraft applications are listed in Table [1] 
These are the Rosenhead-Moore smoothing (Ref. 15]) that yields Scully's swirl 
velocity profile (Ref. (l^), a solid body rotation smoothing (Ref. i29;]) that gives 

^Leishman (Ref. [13|) uses Vg, they are related through Vg = ^^vg. 
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the Rankine swirl velocity profile (Ref. [2^) and a pararaetric Gaussian smooth- 
ing that yields the Lamb-Oseen swirl velocity profile for a = 1.2564312. The 
velocity smoothings for the recently introduced semi-empirical swirl velocity pro- 
file of Ramasamy and Leishman (Ref. [l^l) are also given. The three-dimensional 
velocity smoothing function associated with the Vatistas swirl velocity profile 
(Ref. [27|) for n = 2 is yet to be found. 



2.3 Swirl Velocity Smoothing 

There is another method in use to construct a smoothing from a swirl velocity 
profile. Instead of the two-step process gsip) — >■ 52 (p) — >■ ve{p) one can skip the 
middle step and construct a velocity smoothing directly from a swirl velocity 
profile (Refs. Q and [2^). As an example, Scully's swirl velocity profile can be 
used to construct a three-dimensional velocity smoothing as follows (Ref. j2^ ) , 

2 

-^(^)-;^^53(p) = ^^^(^^. (9) 

If the equivalent two-dimensional velocity smoothing for the right-hand side of 
Eq. [5]is computed using Eqs. |n]and[51 the resulting swirl velocity profile is not 
the intended one, 

arcsinh(-)p 

voip) = , ^ 7 ■ (10) 

A similar exercise can be done for other velocity profiles, but not all of them 
result in closed-form expressions for the final swirl velocity profile. For example, 
when constructing a three-dimensional smoothing from the Vatistas (n = 2) 
swirl velocity profile, 

2 

veip) = ^ g^ip) = —^==, (11) 

the resulting swirl profile is expressed in terms of the generalised hypergeometric 
function 3F2 (Ref. [2]), 

T? (l 1 1-3 5. ^\ 

, ^ 3-^2^2' 2' 4' 4' p4 ^ 

ve[p) = (12) 

which has an extremum of approximately 0.8845 at the non-dimensional distance 
p = 0.7085. 

This method of directly constructing a smoothing from a swirl velocity profile 
is not used in the rest of this paper as it results in incorrect swirl velocity profiles. 



2.4 Velocity of a Viscous Vortex Ring 

When the regularised Biot-Savart law (Eq. ^ is used to compute the induced 
velocity at a point located at the center of the core of a viscous vortex ring, 
the following asymptotic formula in terms of the three-dimensional velocity 
smoothing (73 (p) can be derived (Ref. [11), 



Ua 



log ( — 



47r / g'sip) log pdp 
'0 



(13) 
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with g'^{p) = P^CsCp)- However, this formula is derived under the assumption 
that the shape of the core is circular. More careful analysis shows that the shape 
of the core is a slightly deformed circle which ultimately leads to (Refs. (lOi. i21j) 



r 

AttR 



log ( — 

a 



Jo P Ji P 



(14) 



2.4.1 Examples 



When Eqs.[T3]and[T4]are used with the corresponding three- and two-dimensional 
smoothings, the results are in general not the same, except for the Rosenhead- 
Moore kernel. Substituting its velocity smoothing in Eq. [T3] leads to 



jjRM ^ 



4:TTR 







h(v) 


- 1 







(15) 



where R is the radius of the vortex ring and a the vortex core radius. The same 
result is found when the two-dimensional velocity smoothing is substituted in 
Eq.IH 

For the solid body rotation kernel, the results do not agree. Substituting 
the three-dimensional smoothing into Eq. [13] yields 



^33 



4ttR 



log 



m 

a 



(16) 



Only when the two-dimensional velocity smoothing is used, one gets the correct 
formula as first derived by Kelvin (see e.g. Saffman, Ref. 22|), 



^92 



AnR 



1 I 1 



(17) 



A vortex ring with the two-dimensional parametric Gaussian smoothing has 
a velocity of 



^92 



AnR 



, ,8R\ (I 7 1 2 



(18) 



where 7 is the Euler-Mascheroni constant (Ref. [1]). When the equivalent three- 
dimensional smoothing is used, the result is slightly different, 



r 

4^ 



■o„^)-(:-I.i,„.i 



(19) 



Under the assumption that the core has a circular shape with vorticity dis- 
tribution as given by the Ramasamy-Leishman model, Eq. [T3] can be found to 
give 



R.L 



A-nR 



log(^ 



3 

n=l 



l-5-^log(4M 



(20) 



Values for the constants a„ and 6„ [n £ 1, . . . , 3) are listed in Ref. [18|. The 
formula for the correct induced velocity (using Eq. [T^ contains divergent terms 
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for values of a„ 7^ 1. This means that it is only available for low Reynolds 
numbers, and given by Eq. [TS] where a is replaced with the appropriate value of 
bi. 

Since the three-dimensional velocity smoothing associated with the Vatistas 
swirl velocity profile is not known, no solution is available for Eq. 1131 However, 
the theoretically correct velocity (Eq. [H| for a vortex ring with the Vatistas 
(n — 2) swirl velocity is known, 



AttR 





1" 


h(v) 






~ 2 



(21) 



which is the same as the result obtained with the three-dimensional solid-body 
rotation smoothing (Eg. [T5|. 

In this investigation, it is assumed that the filament core is circular in shape. 
When validating numerically computed induced velocity values, this assumption 
is also used for the reference formulas. Therefore, the formulas derived in this 
section that assume a circular core geometry will be used as reference instead 
of the correct formulas that use no such assumption. The error that is made 
by this assumption can be quantified by computing the difference between the 
resulting formulas for a certain core model. 

For a medium-size transport helicopter in hover, typical values of are 
0(0.1), given that the strength of the tip vortices can be computed with Thov = 
'^'^]'^^f^ (Ref. [il). For the Rankine core, the difference between Eq. [TBI and 
Eq. [T7] is then on the order of 0.025 m/s. A similar difference can be found 
when the two equations for the velocity of a vortex ring with a laminar core are 
compared fEas.fT8land fT9|) . This small difference shows that the use of vortex 
filaments with circular cores in rotorcraft wake simulations does not give a large 
error given that absolute induced velocity values in the wake (in hover) are on 
the order of 11 m/s. 



3 NURBS-based Vortex Filament Method 

Although parametric equations would be sufficient to model the vortex rings 
as used in this paper, a more flexible approach using Non-Uniform Rational 
B-Spline (NURBS) curves (Ref. [13]) was adopted for this research. For el- 
ementary cases such as circles, an exact geometry representation is available 
using rational polynomials. When this is not possible, accurate approximations 



can be constructed easily (see e.g. Ref. [25|). In most filament methods, analyt 



ical solutions of the Biot-Savart law for basic elements of the wake are used. In 
this paper, induced velocity values are computed by numerical integration of the 
regularised Biot-Savart law instead (Eq. [31). Then, adaptive quadrature meth- 
ods can be used to compute values to required accuracy. Or, without changing 
the resolution of the filament geometry, it is possible to trade computing speed 
for accuracy. One disadvantage of the current method is that it can only be used 
when the three-dimensional smoothing associated with a certain swirl velocity 
profile is known (see Table [1]). 
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3.1 Parametric Geometry Representation 

In this section, some background information is presented related to Non-Uniform 
Rational B-Splines. The focus is on concepts that are necessary for a practi- 
cal implementation. In-depth information on the mathematical background of 
NURBS can be found in the book by Piegl and Tiller (Ref. [13 )■ Topics dis- 
cussed here are the constituent parts of NURBS curves: knot vectors and basis 
functions. The equations to calculate point positions and derivative vectors on 
NURBS curves will be shown and finally, a method to construction NURBS 
circles is discussed. 



3.1.1 Knot Vectors 

As generalisations of Bezier curves, (non-rational) B-spline and (rational) NURBS 
curves are constructed by joining together multiple polynomial segments of a 
desired degree to form a single curve. The breakpoints in parametric space along 
the curve where the individual polynomials are joined are called the knots. The 
array of non-decreasing parameters that separate the segments is called the knot 
vector U and has the form U = [uo,ui, . . . , Um] with Ui < Wi+i, i — 0, . . . , to— 1. 

A uniform knot vector has equally-spaced knots. In normalised form, the 
first knot has value and the last knot equals 1. To ensure that a curve inter- 
polates the end points of the control polygon over which it is defined, the knot 
vector is clamped: the multiplicity of the first and last values in the knot vector 
is increased to p + 1, 

U = [ a, ■ -J , a ,Up+i, . . . ,Um-p^i, b, ■ - J ,bj , (22) 
P+i p+i 

where p is the degree of the curve. 

3.1.2 Basis Functions 

The i"^ B-spline basis function of degree p, denoted by Ni^p{u), is defined re- 
cursively by 



1 if Ui < u < Ui+i 
otherwise 



(23) 



Ni,p[u) = N,.p-i[u) H iVj+i.p_i(u). 

Ui+p — Ui Ui+p+l — Ui+l 



When the Biot-Savart law is used to compute the induced velocity at a point due 
to a vortex filament, first derivative vectors are needed at quadrature points on 
the filament. These derivatives can be expressed in terms of the first derivatives 
of the basis functions of the curves, which in turn are defined as 

Nip = iV^p-i(«) — iV,+i,p-i(M). (24) 

Ui+p - Ui Uj+p+i - «,:+! 
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3.1.3 Calculation of Point Locations and First Derivatives on NURBS 
Curves 

A p**^ degree NURBS curve is defined by 

C(w) = (^^N,^piu)w,p}j j ^^7V,,p(M)w;,j a < u < 6, (25) 

with Pi the control points and Wi the weights of the control points. Ni^p{u) 
are the p**^-degree B-spline basis functions defined on the non-periodic and non- 
uniform knot vector U (Eg. E^ . 

Equation can be written as a non-rational, piecewise polynomial curve 
using rational basis functions i?i,p(u), 

n 

C{u)=Y,R^Au)P^ (26) 

where R^Ju) = ^^jj^ff ^^V, ' 

In the non-rational case all weights Wi have unit value and Eq. [51] reduces 

to 

C{u) = N,^p{u)P, a<u<b. (27) 

1=0 

The first derivative at a point u along a polynomial B-spline curve is given by 

n 

C'{u) = Y,N^u)P,, (28) 

1=0 

where the first derivative of the basis function N^p{u) is given by Eq. [511 

Using homogeneous coordinates, NURBS curves can be represented in a com- 
pact and efficient manner as non-rational, piecewise polynomial B-spline curves 
in four-dimensional space using weighted control points P^ = (wiXi, Wii/i, WiZi, Wi), 

n 

C^{u) = J2N,,piu)Pr. (29) 

i=0 

For a (rational) NURBS curve, Eq. [51] can be used directly to calculate the 
derivatives at a point along the curve in homogeneous coordinates 

n 

a-{u)^Y.Kpi^)'PT- (30) 

1=0 

Setting C(u) = ^(")'^('") = ^^r^, where A(u) is the vector- valued function 

^ ^ ^ w{u) w{u) ' ^ ' 

whose coordinates are the first three coordinates of C™(m). Then (Ref [17|, p. 
125), 

A'{u)-w'{u)C{u) 
C [u) = , (31 

w[u) 

where A'(w) is the vector- valued function whose coordinates are the first three 
coordinates of C"^(m) and w'{u) is its last coordinate. 
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3.1.4 Construction of Circles 

Multiple methods can be used to construct circles using NURBS curves (see 
Ref. [Iq). One of the most straighforward approaches uses a set of rational 
quadratic Bezier arcs pieced together by a knot vector with double internal 
knots. As an example, a circle constructed from four rational Bezier segments 
is shown in Fig.[TJ The weights for the even- numbered control points are one and 
the weights for the odd-numbered control points can be found to be Wi = cos 9, 
where 6 is the angle formed by the triplets ZPjPi-iPi+i with Pj one of the 
odd- numbered control points. In case of a four-segment (9-point) circle (Fig.[T|), 
this angle is 45°, so that the weights of the odd-numbered control points equal 
^ w 0.7071. 

3.2 Static Validation 

The numerical method is validated by comparing calculated induced velocity 
values at a point on a vortex ring of unit strength with the analytical results de- 
rived in the previous section. For a range of relative core sizes a/R between 10~^ 
and 1, the induced velocity is computed at the point {—R, 0). The geometry of 
the NURBS-based vortex ring is depicted in Fig. [1] 

The adaptive quadrature method developed by Shampine (Ref. [lij) was 
adapted for induced velocity computations (the original method only works for 
scalar- valued functions, while induced velocity computations are vector- valued) . 
The relative and absolute errors were set to 10~^ and 10"^", respectively. The 
maximum number of subdivision levels is ten and upon subdivision, segments 
are divided in four subsegments. 

Induced velocity values are computed for the core with solid body rotation 
(Eq. [T6| , the Rosenhead-Moore kernel (Eq. [T5|) and lastly, for the Ramasamy- 
Leishman kernel at vortex Reynolds numbers of 10'* and 10^ (Eq. [20]). 

Inspection of Figs. [5] and |3] shows that the numerical and analytical results 
converge towards each other as the relative core size is decreased with a trend 
that is approximately second-order. For relative core sizes smaller than 10~^, 
the second-order trend breaks down. For very small core sizes, the function 
that is integrated starts to resemble finite pulses with a very narrow base. The 
velocity smoothing associated with solid body rotation is the only one that is 
only C° continuous. As a result, the function integrated is not smooth and the 
numerical result breaks down earlier (Fig. l2b| ). These are two examples that are 
difficult for all numerical quadrature methods. 

4 Analysis of Error Sources 

For the segmentation method, two sources of error have been identified in the 
past (Ref. 0). These are errors caused by the incorrect position of straight 
vortex segments (location errors) and errors due to neglecting the local curvature 
contributions to induced velocity values (curvature errors). An additional source 
of error (which was first identified in Ref. [1^) is due to the use of the core 
correction model as introduced by Scully (Ref. [l^). First, the impact of this 
correction on the accuracy and convergence characteristics of induced velocity 
computations will be analysed. Afterwards, the importance of local curvature 
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corrections on induced velocity computations will be investigated. In both cases, 
the NURBS-based vortex filament method is used in the analysis. 



4.1 Influence of Unconditional Use of Perpendicular Dis- 
tance on Accuracy 

The unconditional use of the perpendicular distance in core correction models 
has a detrimental effect on the accuracy of computed velocity of a viscous vortex 
ring. Using the NURBS-based vortex filament method, the induced velocity near 
a vortex ring with the normalised Gaussian velocity smoothing is computed, 
both with the correct smoothing. 



as given in Table [T] and its equivalent that uses the perpendicular distance. 



where p = h/a is the perpendicular distance divided by the core size. 

The velocity induced by the viscous vortex ring is computed in the plane of 
the ring (where only the vertical component is nonzero) on a square with sides 
ten times the ring radius using quadrature rules with eight or sixteen abscissae 
(per quadrant, one rule is used, so induced velocity values are computed with 
either 4 x 8 = 32 or 4 x 16 = 64 evaluations of the Biot-Savart law). Results are 
compared with the analytical solution for a potential vortex ring (Ref. Q). Close 
to the potential vortex ring, the magnitude of the induced velocity increases 
without bounds. When these results are used as a reference for the numerical 
values computed with the viscous ring, it is expected that the errors will be 
large for points close to the ring. For evaluations points in the far field, the 
numerical results should converge to the analytical values. 

The left-hand side of Fig. l4al shows the relative errors w.r.t. the analytical 
results when Eq. [31]is used. The right-hand side of this figure show the relative 
errors when Eq. [33]is used instead. The upper side shows results for m = 8, the 
bottom side shows results for m = 16. A cross section of this data along the 
X-axis is shown in Fig. I4bl Again, the left-hand shows results computed with 
Eq. [32] and the right hand side of the figure shows results for Eq. [33l 

As expected, convergence towards the analytical results is slow near the 
vortex ring. Still, doubling the number of abscissae of the quadrature rule gives 
results with approximately twice as many correct digits. These figures clearly 
show that when Eq. [321 is used as smoothing function, results do not converge 
to the analytical values for evaluation points located on or near tangent lines 
of the quadrature points on the vortex ring. Doubling the number of abscissae 
has the exact opposite effect, for more points in the plane of the ring, results 
fail to converge to the correct values. The width of the regions where results do 
not converge is related to the convergence characteristics of the smoothing used. 
For example, the width for the Rankine solid body rotation model is exactly 
twice the core size, for other smoothings, it is larger. When one would focus on 
comparing results for evaluation points located on the inside of the vortex ring, 
these errors would not show up at all. 




(32) 




(33) 
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When the velocity on the ring itself is computed, similar results are obtained, 
the computed value is reduced with respect to the reference value (Ref. [26{). 
This is illustrated in Fig. [S] where smoothed Biot-Savart function values are 
shown as a function of the azimuth angle of a quadrature point on the vortex 
ring. The area under the curve from to 2tt equals the value of the induced 
velocity at the evaluation point (— i?,0). In this case, the solid body rotation 
velocity smoothing is used with a relative core size of 10%. This particular 
smoothing gives a sharp boundary between smoothed and unsmoothed function 
values. It is clear from this figure that the unconditional use of the perpendic- 
ular distance to correct the Biot-Savart law results in a larger area where the 
potential value is corrected (or reduced) and consequently, gives a value for the 
induced velocity of the ring that is lower than the correct one. 

Using Fig. [SJ the arc length of the smoothing region as a function of the 
relative core size is computed. For the perpendicular smoothing, the width of 
the region is proportional to the square root of twice the relative core size, 

9orig = cos"^(l - (T / R) w y^2a/R > a/R for a < i?. (34) 

When the radial distance is used, the width of the region is proportional to the 
relative core size, 

0new^2sm-\a/{2R))^a/R, (35) 

which confirms the conclusion drawn from Fig. [S] the width of the smooth- 
ing region must be wider for the perpendicular smoothing than for the radial 
smoothing and therefore the computed induced velocity computed with the for- 
mer smoothing is always lower than the one computed with the latter smoothing. 



4.2 Influence of Local Curvature on Induced Velocity 

In methods where curved filaments are discretised using straight segments, the 
two segments directly adjacent to a marker point do not contribute anything to 
the induced velocity at that point. This contribution can be taken into account 
by fitting a circular arc throught the points and using an approximate analytical 



expression for the velocity induced by the arc (e.g. Refs. [23| and [IJI). Using the 
NURBS-based filament method, the influence of a fllament arc at its midpoint 
will be computed for a range of arc lengths and relative core sizes, for various 
smoothings. 

The fllament configuration is shown in Fig. [6] The circular arc with an- 
gle 2A-0 is represented exactly with a quadratic, three-point, rational NURBS 
curve. The control points are labelled Pq, Pi and P2. Induced velocity is eval- 
uated halfway the curve, where it crosses the positive x-axis. Since the curve is 
completely located in the xy-plane, the induced velocity at the evaluation point 
only has a non-zero z-component. 

For smoothings that are at least continuous, results are computed for 
relative core sizes between 0.001 and 0.5 and for filament arc lengths in the 
range 0.1 to 60 degrees. The relative errors between the contribution due to the 
adjacent segment and the velocity due to the whole vortex ring are shown in 
Figs. [7] and m for the Rosenhead-Moore and normalised Gaussian correction. In 
these figures, the relative errors are computed as the absolute difference between 
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the contribution due to the adjacent segment and the velocity due to the whole 
vortex ring divided by the value of the ring velocity, Avz = '^^'^('"'•"''''""^■'•'"g^ 



abs(i)j_^i„g) 

Inspection of Figs. [7] and [5] shows the same general trends, the relative magni- 
tude of the local contribution to the induced velocity at a point on a vortex ring 
is larger for increasingly bigger arc lengths and for smaller core sizes. At dis- 
cretisation resolutions as used in the first generation of wake codes {Aip « 15°), 
neglecting the local contribution results in serious underpredictions of the in- 
duced velocity. For example, an error of approximately 40% is made for a vortex 
ring with a relative core size of 3% for the Rosenhead- Moore smoothing (Fig. [7]) 
when the local curvature contribution is neglected. At a fixed discretisation level 
and relative core size, the relative error becomes smaller for increasing vortex 
Reynolds numbers. For example, at a discretisation of Aijj = 2.5° , the error is 
10 percent for a vortex Reynolds number of one (Fig. [S]). At a vortex Reynolds 
number of 10^, the error is approximately five percent and for a vortex Reynolds 
number of one million, it is even smaller. 

Results shown here were computed with the correct (radial) smoothing, 
not with smoothings that use the perpendicular distance, as most methods 
do (Ref. 0, 23|). It follows that the error contours shown in Figs. [7] and [H] 



cannot be used to get an estimate of the error for segmentation methods that 
use Scully's core correction method. When the perpendicular smoothing is used 
(using for example Eq. [33] instead of Eq. [32]) , errors are orders smaller at fine 
discretisations (compare results in Fig. [S] with those in Fig. [S]). For example, at 
a relative core size of 0.03, the one percent error line shift from Aijj 1° right 
to Alp 10°. From this result, one could draw erronous conclusions regarding 
the need for curvature corrections in a segmentation method. 



5 Convergence Characteristics of the Segmen- 
tation Method 

In this final section, the convergence characteristics of the segmentation method 
will be analysed. The velocity induced by a viscous vortex ring with a certain 
smoothing is computed at three different locations, the first one on the ring 
at (i?, 0), the second one in the centre of the ring and the third one on the 
outside of the ring, at (|i?, 0). In a free wake model, velocities are mainly 
computed at marker points on the vortex filaments. The results presented in 
this section for the point located on the ring itself give a realistic representation 
of the convergence characteristics for these marker points, while the last two 
points are representative for far-field computations. An alternative would be 
to evaluate velocities at a large set of points near the vortex ring for various 
discretisation resolutions, and per discretisation, compute the norm of the error. 
This way unfavourable convergence at a small set of points may be hidden when 
convergence at most points is good. To avoid this, results will be shown per 
point separately. 

The vortex ring is approximated with a set of n straight segments of arc 
length The two modeling options analysed in the previous section are com- 
bined in four separate cases: 1) original core correction without curvature cor- 
rection, 2) original core correction with curvature correction, 3) improved core 
correction without curvature correction and 4) improved core correction with 
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curvature correction. 

The curvature corrections are computed with the NURBS-based filament 
method. Per smoothing, induced velocity values are computed halfway on a 
segment of a circular arc for a range of arc lengths and relative core sizes. By 
storing these values in a two-dimensional array, a simple table lookup can be 
used to fetch the correct value. In more realistic cases, the exact core sizes and 
arc lengths are not known in advance and two-dimensional interpolation should 
be used. 

The original core correction uses the perpendicular distance to an evaluation 
point (Fig. I10a[) to correct the analytical formula for the induced velocity. For 
the improved core correction proposed in Ref. (26l | , the correction factor depends 
on the location of the evaluation point w.r.t. the endpoints of the segment (see 
Fig. llObp . When its projection falls between the end points, the perpendicu- 
lar distance is used as before. If the projection point lies beyond one of the 
end points, the radial distance to the respective end point is used. This way, 
corrections are restricted to a finite area close to the segment. 

For the point located on the vortex ring, the analytical formulas for the 
induced velocity of a viscous ring with a circular core derived earlier are used as 
reference. For the other two points, induced velocity values are compared with 
values computed for a potential vortex ring (Ref. [6|). 

For each smoothing studied, the relative core size a/R was varied between 
0.001 and 0.5 and the number of segments was varied between 5 (A-0 = 72°) 
and 3600 (A'f/; = 0.1°). This covers the complete range of tip vortex core sizes 
and filament discretisations that is used in engineering rotorcraft applications. 
A small, but representative set of these results will be shown and discussed 
hereafter. 

For the point located in the centre of the vortex ring, no curvature corrections 
are applied. In addition, the point is located on the concave side of the vortex 
ring where the use of the improved core correction makes no difference over the 
orginal one. The results for the four different cases all collapse onto a single line 
with a second-order slope. 

Figure [11] shows results for the point on the outside of the vortex ring (at 
(|i?, 0)). Similar to the case where the evaluation point is located on the inside 
of the ring, curvature corrections are zero. Recalling the results from the previ- 
ous section (Fig. Ha]), results may not converge completely to the correct value 
for the original core correction. This is indeed the case: at coarse discretisa- 
tions, second-order convergence is indeed achieved, but for segment arc lengths 
smaller than one degree, convergence stops completely. For the improved core 
correction, convergence is not impeded, even at the finest discretisations. 

Results for the point located on the vortex ring are presented in Fig. [T^ for 
a three percent relative core size, again for the normalised Gaussian smoothing, 
both on an absolute (Fig. I12ap and relative scale (Fig. I12bl) . Second-order 
convergence observed as in Fig. [TT] is nowhere to be seen now. The horizontal 
line at U^^^ « 0.4 m/s is the reference velocity in this case (Eq. [T^ . For the 
four different cases, very different behaviour is observed. 

The original core correction model without curvature corrections signifi- 
cantly underpredicts the velocity at coarse discretisations since it neglects the 
contribution due to the directly adjacent segments. As the discretisation is re- 
fined, the value of the computed result increases until at iS.ip ~ 10° the predicted 
value levels off and converges to approximately 0.24 m/s. As explained before. 
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the original core correction overcorrects induced velocity values which results in 
an underprediction of the total induced velocity. 

When this method is augmented with a curvature correction, results at 
coarse discretisations are significantly improved. The correct value is some- 
what overpredicted due to discretisation errors. As the discretisation is refined, 
the magnitude of the curvature correction is reduced and the errors associated 
with original core correction start to influence the total value. The velocity at 
very fine discretisations converges to the same value as the previous method as 
then they are essentially the same. At the discretisation resolution where the 
computed velocity value crosses the reference value, a small region with a high 
accuracy results. In the plot with relative errors, this shows as a narrow region 
where the relative error has a local minimum. 

For the improved core correction model without curvature correction, results 
at coarse discretisations are exactly the same as for the original core correction. 
As the discretisation is refined, the predicted value keeps increasing until at 
approximately 0.8 degrees, the correct value is crossed. Further refining the 
discretisation results in a further increase of the predicted value, to 0.42 m/s at 
Af/; = 0.1°. From this figure, it is clear that this correction is indeed (just) a cor- 
rection and not correct (see discussion later). However, the error is significantly 
reduced with respect to the original core correction. 

Combining the improved correction with the curvature correction gives the 
best overall results. At coarse discretisations, the curvature correction ensures 
that the predicted value is fairly accurate, while the improved correction pre- 
vents the predicted values from decreasing too much. The curve of the predicted 
values crosses the reference value twice between discretisations of one and two 
degrees. At finer discretisations, the contribution due to local curvature be- 
comes negligible and the values converge to the values as predicted with the 
previous method. 

For a vortex ring with the solid body rotation smoothing and a relative core 
size of 5 percent, results are shown in Fig. [T21 In this case, the smoothing is 
not continuous which results in nonsmooth covergence characteristics. The 
positions of the kinks in the curves (both for the original and the improved 
core correction model) can be predicted analytically due to the sharp boundary 
between corrected and uncorrected regions for this smoothing. 

For the original core correction model, the following implicit equation gives 
the relationship between the relative core size and the arc length of the segments 
that results in the position of the kinks, 

At/. (2n + l)A^ a 
cos— -cos = -, (36) 

where the parameter n is the n"^ adjacent segment to an evaluation point on 
a filament that has a nonzero contribution at that evaluation point. Eq. 1361 is 
implicit in Atp and must be solved iteratively. A good initial guess for Aip = 
2^j^i arccos (l — ■^). As an example, for n = 1 and n = 2 (which correspond 
to the first and second adjacent segment with a nonzero contribution at the 
evaluation point), the values for Aij; equal 12.88°and 7.43°, respectively. 

For the improved core correction model, the equation for the relationship 
between the relative core size and the arc length of the segments that gives the 
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position of the kinks is explicit, 

AV' = - arcsin (^) . (37) 
n \2R/ 

For the first two points, this gives values for A^/; of 2.87°and 1.43°. Tick marks 
are added in Fig.[T3]at these points. 



6 Discussion 



The correction proposed in Ref. 26| was shown to give improved convergence of 
the segmentation method for resolutions to approximately one degree segments, 
but it is still just a correction. When even finer segment discretisations are 
desired, the exact correction may be needed. It can be found by correcting the 
Biot-Savart law with a velocity smoothing before integration instead of correct- 
ing it afterwards with the desired swirl velocity profile. As an example, ana- 
lytic integration of the Biot-Savart law regularised with the Rosenhead-Moore 
smoothing for a straight line segment, 

RM _ r /■ r dr{u) 



will give the exact correction for Scully's swirl velocity profile. This should then 
be repeated for each smoothing one wishes to use. Using the NURBS-based 
filament method, it is easy to show that if one would use the analytic result of 
Eq. [3S] with a segmentation method, second-order convergence can be achieved 
for points on the vortex ring. An exact NURBS representation of a segmented 
vortex ring (a regular polygon) with n segments is created and the induced 
velocity is computed at a point on the ring using the regularised Biot-Savart 
law. Results are shown in Fig. [T3] for the normalised Gaussian smoothing with a 
three percent core size both on an absolute and relative scale, for both the case 
without curvature correction and the case including curvature correction. Re- 
sults for the segmentation method using the improved core correction similar to 
the results shown before (Fig. [T2|) are included for easy comparison. Inspection 



of Fig. Il4bl shows that for the case without curvature correction at discretisa- 
tions finer than Ai/; w 3°, a second-order trend is observed. The addition of 
a curvature correction improves the convergence of segmentation methods at 
coarser discretisations. Results computed using both methods converge to the 
same values for fine discretisations. 

Derivation of the correct, smoothed version of the velocity induced by a 
straight line segment will only improve the accuracy of results at discretisations 
finer than what is currently used in engineering wake models, so it may not be 
worth the effort. Furthermore, the derivation may become more difficult when 
the vortex strength or the core size is allowed to change along the filament. 

Bhagwat and Leishman (Ref. Q) have computed induced velocity values at 
points near and on a potential and viscous vortex ring for a very fine discretisa- 
tion and used that as a reference value to compute convergence characteristics 
at coarser discretisations. The conclusion of Bhagwat and Leishman is that: 

"Therefore, the discrete straight-line approximation for the in- 
duced velocity of either a potential or viscous vortex ring is formally 
second-order accurate." 
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However, one of the assumptions used to arrive at this conclusion, that the 
result at the finest discretisation is the correct one and that it can be used as 
a reference value to compute errors at coarser discretisations, is not true. As 
shown here, the straight-line approximation method can only attain second- 
order convergence for points on the ring provided that formulas are corrected 
before integration, and not afterwards. 

7 Conclusions 

The effect of two modeling assumptions used in straight-line segmentation meth- 
ods have been examined rigorously. The first one is related to the assumption 
that contributions due to local curvature are sufficiently small so that they may 
be neglected completely. The second assumption is related to the method used 
to regularise the analytical solution for the induced velocity due to a single, 
straight vortex line segment. 

A high-order vortex filament method is presented to investigate these model- 
ing assumptions. The method uses Non-Uniform, Rational B-Spline (NURBS) 
curves as a basis for filament geometry representation. The Biot-Savart law is 
regularised with a velocity smoothing function and induced velocity values at 
evaluation points are computed using an adaptive quadrature method. The re- 
lationship between three-dimensional velocity smoothing functions and the asso- 
ciated two-dimensional swirl velocity profiles was shown. The three-dimensional 
smoothings related to various swirl velocity profiles were given except for the 
swirl velocity of Vatistas (n = 2), for which it has not been found yet. The 
method is validated statically by computing the velocity of a viscous vortex 
ring. For decreasing relative core sizes, numerical results arc shown to converge 
to the value predicted with approximate analytic expressions that assume a cir- 
cular core shape. These are different from the (correct) expressions where no 
assumption is made regarding the core shape, except for the Rosenhead-Moore 
velocity smoothing. 

The NURBS-based vortex filament method is subsequently used to study the 
effects of two common modelling assumptions used in straight-line segmentation 
methods, the use of core corrections and curvature corrections. For the former 
one, the original model and an improved version are tested. This analysis shows 
that for evaluation points located near tangent lines of quadrature points, the 
original core correction model predicts induced velocity values that are too low. 
For points on the vortex ring, this results in a consistent underprediction of the 
ring velocity by approximately 40%. 

The consequences of neglecting the contribution from the adjacent part of a 
curved filament to the induced velocity are also studied using the NURBS-based 
vortex filament method. It is shown that for a range of segment discretisa- 
tion resolutions as used in engineering rotorcraft applications, discarding these 
corrections gives errors between 1 and 40 percent, depending on the relative 
core size, the arc length, the particular smoothing and optionally, the vortex 
Reynolds number used. When the original core correction is used, these errors 
are orders of magnitude smaller for the same relative core size and arc length. 

Finally, the segmentation method is analysed in detail. Induced velocity is 
computed at three points in the vicinity of a viscous vortex ring. The first one 
is located at the centre of the ring, the second one is located on the outside 
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and the last point is located on the ring. Results for the former two points are 
representative for far-field computations while results for the latter point are 
representative for self-induced velocity computations. It is concluded that for 
points in the far-field of the vortex ring, second-order convergence can only be 
achieved with the improved core correction model. For points on the vortex 
ring, the best possible results over a wide range of discretisations are achieved 
by combining the curvature correction with the improved core correction. The 
former ensures that computed values are close to the reference for coarse dis- 
cretisations while the latter gives good results at fine discretisations. For the 
improved core correction at discretisations finer than 5 degrees, the use of cur- 
vature corrections gives results that are an order of magnitude more accurate. 
When the discretisations are finer than one degree, differences are so small that 
curvature corrections may be neglected. Finally, it is demonstrated that the 
use of corrections applied to the Biot-Savart law after analytical integration 
precludes second-order convergence to the correct value. This is only possible 
when the Biot-Savart law is regularised before integrating it analytically. 
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P3 P2 Pi 




Figure 1: Geometry of the NURBS circle as used to compute vortex ring induced 
velocities. 
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Table 1: Three- and two-dimensional velocity smoothings and swirl velocity profiles relevant to rotorcraft appli- 
cations. 



Smoothing model 53 (p) g2{p) ve{p) 

Rosenhead-Moore^ 2! — 



Solid Body Rotation, p < 1 
Solid Body Rotation, p > 1 

Parametric Gaussian^ 
Vatistas {n=2f 



4.(p2 + l)f ^ 



a(p)-p-\/l-p^ 



271-2 27r 
J_ J_ i 

eri{p^/a)-2p^ exp(-ap2) l-e^p{-ap'^) l-exp(-ap2) 
47r 27r p 



1 an(erf(p\/b;r)-2pV^exp(-b„p2)) 1-En=i °" exp(-6np2) 1-En=i g^P(-'''»P^) 



Ramasamy-Leishman ^"=^ ""^^^ ~ ^^"='1; 



^ Equivalent of Scully swirl velociy profile. 

Equivalent of Lamb-Osecn swirl velociy profile when a 1.2564312. 
^ Three-dimensional smoothing not known. 



P 





^0-4 lQ-3 lQ-2 

Relative core size a/R 
(b) 



Figure 2: Viscous vortex ring velocity as a function of relative core size a/R for 
the smoothings of Rosenhead- Moore and Rankine, (a) absolute values and |(b)| 
relative errors. 
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Figure 3: Viscous vortex ring velocity as a function of relative core size a/R 
for the smoothing based on the Ramasamy-Leishman swirl velocity profile at 



vortex Reynolds numbers of 10 and 10 , (a) absolute values and |(b)| relative 
errors. 
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Figure 4: Convergence of induced velocity values computed in the plane of a 
viscous vortex ring with the normalised Gaussian smoothing using the NURBS- 
based vortex filament method. Results on the right-hand side are computed 
with the perpendicular core smoothing, while results on the left-hand side are 
computed with the correct, radial smoothing. In (a) results for quadrature 
rules with 8 abscissae are shown in the top half while those for 16 abscissae are 
shown in the bottom half. A cross section of the data along the x-axis is shown 
in[(b)l 
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Figure 5: Boundaries between potential and corrected (viscous) Biot-Savart 
function values for a point {—R, 0) on a vortex ring {a/R — 0.1) for the Rankine 
swirl velocity profile. For both corrections, the area under the curve from 0° to 
360° in (a) equals the induced velocity at the evaluation point due to the whole 
vortex ring. The width of the correction regions are 9orig — 25.8° (using Eq. 



and 6ri 



5.7° (using Eq. 
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Figure 6: Configuration of the portion of the vortex ring that is not taken into 
account when computing induced velocities at point P. The induced velocity is 
computed from the solid part of the circle by substracting the local contribution 
from the velocity of the complete vortex ring. 




Filanicnt arc length 2Aip [deg] 



Figure 7: Relative error in predicted velocity at a point on a vortex ring with 
the Rosenhead-Moore velocity smoothing by neglecting the influence of the con- 
tribution of the adjacent part of the ring. 
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Figure 8: Relative error in predicted velocity at a point on a vortex ring with 
the normalised Gaussian velocity smoothing by neglecting the influence of the 
contribution of the adjacent part of the ring. 
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Figure 9: Relative error in predicted velocity at a point on a vortex ring with 
the normalised Gaussian velocity smoothing by neglecting the influence of the 
contribution of the adjacent part of the ring for the perpendicular smoothing, 
using Eq. [331 
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(a) Original core correction. 




(b) Improved core correction. 



Figure 10: Two core corrections that can be used to correct the velocity field 
induced by a finite, straight vortex line segment. The gray region is the area 
where the correction is applied. Its width depends on the convergence charac- 
teristics of the specific smoothing, for the solid body rotation smoothing, c = 1, 
for other smoothings, c 1. 
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Figure 11: Convergence of straight-line segmentation method at (0, 1.2R) for a 
viscous vortex ring with Gaussian smoothing {a/R = 0.03). Convergence is not 
satisfactory for the original core correction model. 
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Figure 12: Induced velocity computed with the segmentation method for a point 
located on a viscous vortex ring with Gaussian smoothing {a/R = 0.03): (a) 
absolute values and |(b)| relative errors. 
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Figure 13: Induced velocity computed with the segmentation method for a point 



located on a viscous vortex ring with Rankine smoothing (cr/i? = 0.05): (a) 
absolute values and |(b)| relative errors. Positions of the labels at the upper side 
are computed with Eg s. [Ml and 1571 
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Figure 14: Induced velocity computed with the segmentation method for a 
point located on a viscous vortex ring with Gaussian smoothing (cr/i? = 0.03) 
using the improved core correction model (segmented ring) and the numerically 



computed exact correction (using a NURBS-based regular polygon) : (a) absolute 
values and 1(b) I relative errors. 
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